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Abstract 

Investigations are made on the saddle point calculations (SPC) under the auxiliary 
field method in path integrations. Two different ways of SPC are considered, Method(I) 
and Method (II), to be checked in an integral representation of the Gamma function, 
r{N), as a bosonic example and in a four-fermi type of Grassmann integral where one 
"fermion mass" ujq differs from the other iV-degenerate species. The recipe of Method(I) 
seems rather complicated than that of (II) superficially, but the case turns out to be 
opposite in the actual situation. A general formalism allows us to calculate for r{N) 
up to 0(l/iV^^). It is found that both happen to coincide in the bosonic case but in 
the fermionic case Method (II) shows a huge deviation in the weak coupling region where 
Wo <C 1. 
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1 Introduction 



The auxiliary field method (AFM) is one of the most powerful approximation scheme in 
path integrations. The recipe for AFM is given as follows: suppose a partition function, 



Z = / d^ad^a* exp 



A2 2' 



N 



(cr* • cr) = ^ CJ*CJa , (1) 



where cJa's, (a = 1, . . . ,N), are fermionic or bosonic degrees, u) (called a mass) and (a 
coupling constant) are parameters. 

(i) Introduce the auxiliary fields [H |2], y , (or Hubbard-Stratonovich Field in the solid 
state physics[3]) by inserting the identity in terms of the Gaussian integral, 



Z= d^ad^a* e— (-*-)-^(-*-)' 



dy 



X / exp 



in order to remove the four-body interaction, to obtain 



A 



z 

(ii) Put y ^ \/]V y and perform the "Gaussian" type integration of da to write 
E{y) = j duadal exp [- (u; + zAy) alaa] ; for . 



Thus Z oc i?^, yielding 



, (2) 



(3) 



(4) 



(iii) Write 



-N(^-lnEiy) 



(5) 



In E{y) , 



(6) 



and assume 1— > cx) to perform a saddle point calculation (SPG): find saddle point(s) yo, 
satisfying the stability condition. 



f'iy) 



, /'•^■'(yo) > ; (stabihty condition) , 



(7) 



and expand f{y) around yo, which gives us a power series of called the loop 

expansion[3]. This is not a convergent but an asymptotic series, of course. 

The prescription is simple and straightforward compared to other nonperturbative 
methods such as the variational [5] and the optimized perturbation [6]. Moreover, the 
studies in 1- dimensional(= a quantum mechanical) as well as 0-dimensional(= an inte- 
gration) bosonic[7] and fermionic [SI |9] models tell us that we can obtain a fairly accurate 
result, even w^hen A^ = 1 or small, from the weak to the strong coupling A, by taking 
higher loops into consideration properly. 
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Now consider a slightly generalized model; 

a ad (7 exp 
where {N + 1) x (A^ + 1) matrix u is given 



icr* ■ u ■ cr) (a* ■ cr) 



(8) 



N 



with In and being the N x N unit matrix and the dimensional zero vector. This 
corresponds to a field theoretical model of interacting (Tq and da (a = 1, 2, . . . N) bosons 
or fermions with the "mass" ujq and uj respectively. Introduce an auxiliary field, y, as the 
above, to find 

= j^dyg{y) exp [-Nf{y)] , (10) 



where 



9{y) = j daedal exp [- (wq + «Ay) crQcro] , (11) 

and f{y) has been given by ([6]), with 

E{y) = j duadal exp [- (w + iAy) o-*crJ ; for . (12) 

In view of ([TO]) an issue comes up: there are two ways of performing SPC. 
• Method(I): iifl, 



dt5(i)e-^^W , (13) 



find the saddle point of /(t), /'(to) = Oi then expand f{t) as well as g{t) around 
to- Here and hereafter we adopt t instead of y as the integration variable. 

• Method(II): rewrite ^ as 

In^ j (ite-^/W , fit) ^ fit) - 1 Ingit) , (14) 

then find the saddle tc of fit), f'itc) = 0, and expand fit) around tc- Finally put 
tc, given in terms of 1/A^ series, into the expression. 

If becomes large both results would match but, as mentioned above, our interest 
is to study the validity of AFM when A^ is small. In this paper, therefore, we study 
the difference between two methods by considering bosonic and fermionic integrations 
(0-dimensional field theoretical models). In SeclH we develop a general formalism of SPC 
and calculate the asymptotic expansion of the Gamma function as a bosonic model in 
Seci3l In SecSl we examine a four-fermi type Grassmann integral and the final Secl5] is 
devoted to a discussion. 



^The integration range of t, in (fT5|) and (|14p . need not be specified. Those are given from — oo to oo in the 
final expression under the saddle point method. (See (|15p to ^TE\i .) 



3 



2 Saddle Points and the Asymptotic Expansion 

In this section, we develop a general formalism of SPC. 

• Method(I): start with the expression (|13p and expand all the integrands around the 
saddle point, Iq, to find 



Im = e-^^« 



oo (m) 



/^ME^(*-*or exp 

\m=0 ■ / 



^(2) 

{t - tof - Rn 



(15) 



where 



JO 



n! 



n=3 



with 



(16) 



(17) 



After checking the stability condition ([7]), /q^^ > 0, put t— to = x/y NfQ^\ while assuming 
N large in p^ . to find 



-Nfo 



IN 



E 



(-)^5o ^ 



m/2 



dxx"* (i?7v(x)) e 



A: „-a;2/2 



(18) 



where (from (fT6 



-j^ oo 



Fn+3 



JO 



in) 



A^(n + 3)!iV"/2 ' 

n=0 



(2) 



n/2 



(19) 



(The integration range of x now stretches from — oo to oo.) Therefore 

(m) • . i'™ I I. I T^fc 



/TV 



/~7(^ ^ m\k\ / ,(2)\"V2 



1 ^^(m+fc+E,tl'^J)/2 

^ i A 



^ni+3 /n2+3 



rifc+S 



,n2,„rn.=0 
oo 



Now put 



(ni+3)!(n2 + 3)! (nfc + 3)! 

m + A; + X]j=i "-i 



(20) 



L = 



f211 



which must be integers, L G Z, not half- integers, L G Z + 1/2; since the power of x in 
(1201) reads 



+ 3/c + ^nj = 2(L + /c) , 



(22) 
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so that the integral vanishes unless L -\- k £ leaving us 

/oo 
dx e-"'/2 = ^(2(L + k)-l) 
-oo 



The expression ()20p then turns out to be 

oo , 2L 



27T 1 ^(-)^ 



0^ 



where 



F(ni,n2, . . . ,nfc) 



L=0 fc=0 

X ^ F{ni,n2, . . . ,nfc) , 

all possible {nj} 



(21 -k- Ej=i /^(2)^ (2i"fc-E,ti n,)/2 U (n, + 3)! 



(23) 



(24) 



, (25) 



is a symmetric function of nj's. In ()24p the sum should be taken under the condition, 



Uj < 2L — k ; 

since m > in ([22]). Note if /c = 0, then X]°=i = and n°=i G{nj) = 1, so that 



(26) 



F{ni,--- ,nk) 



k=o 1 go 



(2L) 



(2L)! /^(2)^^ ■ 



(27) 



The conditional sum of n^-'s, (j26p . can be expressed as an alternative form: suppose Qc 



of rij's are alike of A^, Uj-^ = nj^ 



rijQ^ = Aa, to write 



/ Qi 



F 



Q2 



Qp \ 



.(28) 



It is clear that the condition (l26l) reads as 



P<k ; ^Q^ = k, ^QaAo,[=^nA <2L-k . 



(29) 



a=l 



The multiplicity reads 



k \ ( k — Qi 



Qi J \ Q 



a=l 



Qp 



k\ 



(30) 
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so that becomes 

2L 



V Jo L=0 k=0 



j{2) 

X E 5^<K}.K"}.-.K}) . (31) 

all possible {Aa} 
ELi Q-=fc ; (P<fc) 

where we have assumed that Ai < A2 < ■ ■ ■ < Ap. This is the main formula of the 
Method(I). Let us classify 1^, according to the WKB-approximation [1] under which 1/h 
appears instead of N; 

1. Tree: 

(/7v)tree = e-^/°<70 ■ (32) 

2. Woop: terms up to / (> 1) in (pi|) : 



(^-)Moop - e--/" ^^,0 (l + • • • + O (^) ) . (33) 

• Method(II): by putting f[t) i— > f{t) ; (^(t) i— > 1 and to ^ ^c; all expressions in Method(I) 
can be read as those of Method(II). Write 

/W^/W(t,) , (34) 

then (1211) is changed to 



2L k Ej*=i"j=2i-fc 

^E^(2(L + A:)-1)!! F{nun2,...,nk) , (35) 

'==0 all possible {rij} 

k p ^(n) 

^^=i(^i + 3)! f/P) 

It should be noted that without g{t) the conditional sum is given by the equality, 
X]j=i '^i = 2L — A;, not by the inequality (|26p . Also (j3ip is changed to 



^ E o^^(K*}.K}.-.K}) . (3^) 

all possible {Aq} 
ELi Q«='= ; {'P<fc) 

where , j^f '| , . . . , j^p^}) is defined by the expression, ([28]), with the tildes. 

Here again the conditional sum (j29p becomes simpler, ^^^j^ Qa^a = 2L—k. Classification 
in this case is again 
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1. Tree: 



(/;v)t.ee = e-^/^-' , (38) 



(0) 



2. Moop: terms up to I (> 1) in ([20]) or i^: 



(Whoop -e-/?y-^(l + ... + o(^)) . (39) 



Note that there are additional powers of hidden in tc- The final task is then to 
expand all functions of tc up to (1/A^)'~^, which however depends on individual 
models so is relegated to the following sections. 



3 A Bosonic Case: the Gamma Function 

As a simple bosonic example, in this section we consider the Gamma function, 

dtt^-ig-* = N^In ; In = J dt (^ij e-^^*^'-^*) , (40) 

where we have put t Nt in the final expression. Then 

Method(I) : f{t) = t-\nt- g{t) = ^ . (41) 
Method(II) : f{t) = t - M - j Int . (42) 



»Method(I): the saddle point, to is given by 

1 



f'{to) = 1 - 1 = ^ to = 1 , (43) 



so that 

/o = 1 = 50 , ft^'^ = -gt^ = -i-rm ; (n > 1) , (44) 

(2) 

which ensures the stability condition, /q = 1 > 0, and 

k ^ 

F(ni,n2,...,nfc) = TT- — — , (45) 

from ([25]). Therefore from ^ and ([3T]) . 

rz oo ^ 2L 

L=0 fe=0 
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where 



T{L,k\<2L-k) ^ i^(2(L + A:)-l)!! 11 (^TTs) ^'^^ 

all possible {nj} J-^ 



(_)^(2(L + A:)-1)!! 



Qi!---gp! 

all possible {Ac} 



ELiQc=fc ; iP<k) 
1 



(^1 + 3)^^ {A2 + ■■■{Ap + 3)^^ 



(48) 



Here we concentrate on the latter expression (j48p and realize that most of the terms in the 
conditional sum, X]r=i Qa^a < 2L — k, are canceled, leaving only the terms satisfying 
the equality J2a=i Qa^a = 2L — k , that is 



2L 2L 



^ r(L, k\<2L-k) = ^ T{L, k\2L - k) , (49) 



fc=0 fc=0 



to give 



In 



0^ ^ , 2L 

i-~''Y.wY.nL.m-k) , (50) 

L=0 A:=0 



where 



r(L,A.|2L-fc)^i^(2(L + A:)-l)!! 11 (^1) 

all possible {uj] 

= (-)'= (2(L + ^^)-l)!! 



Qi\---Qpl 

all possible {a^} 

Er=i Qc=fc ; (P<fe) 

X (52) 

(Ai + 3)^1 (^2 + 3)^2 • • • (^P + 3)^^ ' 

with the definition when k = 0, (implying P = 0), and L 7^ 

T(L,0|2L)=0, (53) 

since otherwise there is inconsistency in the conditional sum X]r=i Qa^ai= 0) and 2L{^ 
0). The numerical values of T{L, k\2L — k) (apart from the factor (2(L + fc) — 1)!!) are 
listed in the appendixEl The proof of ([M]) is rather lengthy then relegated to the appendix 

m 

In view of (j24p and (I25p . the above fact implies that there is no contribution from 
g^\m > 1), that is, the equality holds 



00 

dt 
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under the expansion (Actually it does hold without 1/N in this case: see the expres- 
sion ()167p . A detailed discussion is relegated to the appendix iBl) 
In view of ^ as well as (|iO|) . 



(^(^))trcc 

And the result of 15-loop is from ()33p . 



(55) 



(r(iv)) 



15— loop 




2 L 

^ 12N^ 288 A^2 



139 1 



571 



1 



5246819 



1 



51840 iV3 2488320 iV^ 
534703531 1 



209018880 
4483131259 



+ 



86684309913600 
6232523202521089 



75246796800 902961561600 iV^ 

1 432261921612371 1 

+ 



514904800886784000 A^^ 
1 25834629665134204969 



1 



86504006548979712000 
1579029138854919086429 



13494625021640835072000 A^^^ 
1 746590869962651602203151 



+ 



1 



+ 



9716130015581401251840000 N^^ 116593560186976815022080000 N^^ 
1511513601028097903631961 1 



2798245444487443560529920000 N^'^ 



+ 0(N' 



(56) 



This coincides the results in referencefl; up to 1/N^^U\, and 




N. 1 - 
N= 2 



Figure 1: The ratio of Approximation to Exact in Method(I): the horizontal axis designates L, 
< L < 14, omitting the tree part. The solid line with the circle, the dotted with the triangle, 
the dashed with the cross, and the dash-dotted with the square designate = 1, 2, 5 and 10 
respectively. The tail end of the A^ = 1 line deviates from the unity, implying the asymptotic 
character of the 1/A^ expansion. 



In the table [H we list A^ = 1, 2, 5, and 10 results up to 15-loop including the tree ones, 
where the optimized values are shaded. Also we plot the ratio of approximate to exact 
values for < L < 14 (omitting the tree part) in the figure [TJ From these we convince 
the validity of the loop expansion in this case; since even in the smallest A^ = 1 case, 
the 2- or 3-loop approximation gives ~ 0.1 or ~ 0.2% error. It also should be noted that 
the characteristic feature of the asymptotic expansion can be read from deviation after 
passing through the optimized values, which is most clearly seen in the figure [T] at the tail 
end of the A^ = 1 line. 

^We have no information how those results are obtained, but the derivations would be from another pre- 
scription such as the Watson's Lemma [13] 



9 





N" — 1 N — 2 N — b N — 10 


Exact 


1 1 24 362880 


tree 
(tree) /Exact 


0.36788 0.54134 2.10561 X lO"^ 4.53999 x 10^ 
0.36788 0.54134 0.87734 1.25110 


1-loop 

(1-loop) /Exact 


0.92214 0.95950 2.36038 X 10^ 3.59870 X lO'^ 
0.92214 0.95950 0.98349 0.99170 


2-loop 
(2-loop)/Exact 


0.99898 0.99948 2.39972 x IQ-"- 3.62868 x 10^ 
0.99898 0.99948 0.99988 0.99997 


3-loop 
(3-loop) /Exact 


1.00218 1.00031 2.40005 X 10^ 3.62881 X 10^ 
1.00218 1.00031 1.00002 1.00000 


4-loop 
(4-loop) /Exact 


0.99971 0.99999 2.40000 x 10^ 3.62SS0 x 10^ 
0.99971 0.99999 1.00000 1.00000 


5-loop 
(5-loop) /Exact 


0.99950 0.99998 2.40000 x 10^ 3.62880 x 10^ 
0.99950 0.99998 1.00000 1.00000 


6-loop 
(6-loop) /Exact 


1.00022 1.00000 2.40000 X 10^ 3.62880 x 10^ 
1 .00022 1 .00000 1 . 00000 1 . 00000 


7-loop 

(7-loop)/Exact 


1.00029 1.00000 2.40000 X 10^ 3.62880 X 10^ 
1 .00029 1 .00000 1 . 00000 1 . 00000 


8-loop 
(8-loop) /Exact 


0.99974 1.00000 2.40000 x 10^ 3.62880 x 10^ 
.99974 1 .00000 1 . 00000 1 . 00000 


9-loop 

(9-loop) /Exact 


0.99969 1.00000 2.40000 x 10^ 3.62880 x 10^ 
0.99969 1.00000 1.00000 1.00000 


10-loop 
(10-loop)/Exact 


1 00047 1 00000 9 Aonnn v in-'- fiossn v 
1.00047 1.00000 1.00000 1.00000 


11-loop 
(11-loop) /Exact 


1.00053 1.00000 2.40000 X 10^ 3.62880 X 10^ 
1.00053 1.00000 1.00000 1.00000 


12- loop 
(12-loop)/Exact 

13- loop 
(13-loop)/Exact 


0.99877 T.OOOOb 2.40000 X 10^ S.'eis's'O x 10^ 
0.99877 1.00000 1.00000 1.00000 
0.99862 1.00000 2.40000 X lO'^ 3.62880 X lO" 
0.99862 1.00000 1.00000 1.00000 


14-loop 
(14-loop) /Exact 


1.00452 1.00000 2.40000 X 10^ 3.62880 X 10^ 
1.00452 1.00000 1.00000 1.00000 


15-loop 
(15-loop) /Exact 


1.00502 1.00000 2.40000 x 10^ 3.62880 X 10^ 
1.00502 1.00000 1.00000 1.00000 



Table 1: Results of Metliod(I) in iV = 1,2,5 and 10 up to L = 14. Even in iV = 1, 2- 
loop(L = 1) approximation is sufficiently close to the exact value, however, deviation becomes 
gradually eminent in the loops larger than 7. The optimized values are shaded in each N, 
whose position shows that the 1/A^ expansion is indeed an asymptotic one. 



•Method(II): from (H2]) . the saddle point is 

fit) = 1 - (i - ^) 7 = ^ = ^0 + 77 ; to = i,ti = -i 

and satisfies the stability condition, 



U N-l 



(57) 



(58) 



Note that 



to give Fn 



/, = ie(l-lnt,) , /W^/»(t,) = (-)-(n-l)!— 1^ ; (n > 2) , 



n/2-l 



SO that 



(tc 



(59) 



, (60) 



with rij + k = 2L, from (f36l) . Therefore with the aid of (f35]) 

2L 



In ~ e" 



~Ntc{l~lntc) 



27rtc 

L=0 



k=0 



(61) 
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The tree and the 1-loop part are given 



{I. 



-Ntcil-lntc) 



-N 



ti=0 



(-^iv)i_loop 



-Ntc{l-lntc) I'^^'tc 



N 



-N 



ti=0 




(62) 
(63) 



From the 2-loop approximation, ti([57D should be included to tc, and all function of tc 
must be expanded up to 0{1/N^) : L = 1,2 . . . : in (j6ip . introduce the prefactor V, 



V 



-Ntc.{l-\ntc) 



2'Ktc 



-N 




1 1 



1 1 



12 iV2 12 ATS 1440 Ar4 + I ^5 



,(64) 



and the loop factor C 



L=0 ' k=0 

7i{l) + 2n(2) + n(3) n(l) + 3n(2) + 3n(3) + n(4) 



+ 



ATS 



+ 



where 



2L 



n 



(L)E^ J]r(L,A;|2L-A:) ; 

k=0 



n(l) = 1 , n(2) 



1 



n(3) 



139 



, n(4) 



571 



288 ' ' 51840 ' ' 2488320 
from the table [5] in the appendix [Xj V x C gives the 5- loop approximation of In 

(^^)5-loop 



(65) 

(66) 
(67) 



-N 




L=l 



which is exactly the same to the one (L i-^ 4 in ([50|) ) in Method (I). Many terms in 
the numerator, n(l),...,n(L — 1) are canceled, leaving us only n{L). These miracle 
cancellations occur for all orders of 1/A^, yielding the result that there is no difference 
between Method(I) and (II) in this case. 

The reason is rather easily figured out: because of the formula (|167p in the appendix [Bl 
we can put g{t) = 1/t in the integral In (j40p to the unity, g{t) i— > 1. In other words, 
Method(I) is equivalent to Method(II) in the Gamma function case. 



4 A Fermionic Case 

Our target is a fermionic version of dH); 

Z = J d^'+^id^'+^i* exp 
where ^,^*(^,^*) are + l(A^)-dimensional Grassmann variables. 



(69) 



(70) 
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with 



^7V+i^" ^ d^^d^o ^ d^Nd^N^i ■■■dCo, d^'+'e ^ dCod^'e ^ dCodCi • • • , (71) 

and (A^ + 1) x (N + 1) matrix, u, has been given in ([9]). Z is calculable by means of a 
standard Grassmann integration; 



/ 



' n\5r, 



(72) 



to obtain 



[fi 



2 = E 



r=0 



(73) 



In this analysis, we assume that aU parameters in this model are real and positivH, 
w > 0, Wo > 0, A > 0, and take N = 2 with 



< A < 10 ; uJo = W^LO , uj , lO^^u; . 

(The case, loq = lO^to is a toy model of u-, d-, and s-quarks.) 
Introducing an auxiliary field, in terms of 



dy 



1 

-oo \/2vr 

into the target (f69]) . we obtain 



exp 



Z = J — dy{ujQ + \y) exp 
V 2vr 

where we have performed the Grassmann Gaussian integration. 



and y has been scaled, y i— > v Ny, as before. Now write 



Z 




'N ; 



In = j dt g{t) e" 



■Nf{t) 



for Method(I) and 



g{t) = wo + At ; f{t) 



In= dt e-^^'W 



In (w + At) , 



ht) = f{t) - ^ \ng{t) = 1 - In (w + At) - 1 In {u, + At) 
for Method(II). (Here y has been switched to t.) 



(74) 



(75) 



(76) 



(77) 



(78) 
(79) 



(80) 



^Although we can see an interesting phenomenon when A^ < ;the caustics emerge[9]- 
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Method(I): the saddle points are determined by 



t=to U> + Ato 

Here and hereafter the equation is called as the gap equation |14j. If we introduce 

no = co + Xto , (82) 

(jHT]) becomes 

{Qof - ujHo - = , (83) 

yielding to two saddle points 

The stability condition ([7]) in this case reads 



which is positive for both ^^q^"*- The value of f{t),g{t) and derivatives at the saddle points 
are given by 

fo=^-^^^^-lnno ; ft^ = {m-iy.^-^y ; {m>3), (86) 
go = f2o + 5LO ; = ^ ; 5o"^ = ! i^n>2) , (87) 

where 

6uJ = UJQ — UJ 

Then from ()19p 

-^n,+3 _ 1 ^^{2)\-K+3)/2^(n,+3) _ 1 / -g(^o) A 



where 



is the sign function. Accordingly, in view of ()24p with (j25p . we find 

^(2(L + A:)-1)!! ^ F(ni,n2,...,nfc) 

all possible {nj} 

2(L+A:) 

[(r2o + 5'^)r(L,fc|2L-A:)-r?oT(L,A:|2L-fc-l)] , (91) 



(W + A^ 
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where use has been made of the notation (j5ip : T{L, k\2L — k — 1) is defined by replacing 
the sum X]j=i = 2L — /c to X]j=i 'b' = 2L — A; — 1 in (I5ip . Therefore from (|25p 

2(L+fc) 



: exp 



-iV 



2A2 



A 



oo ^ 2L 

A/T , 
L=o'^ fc=0 \ A/(f?o) + A2 



X [(I?o + <^^) T^l^, ^|2L - A:) - noT{L, k\2L - k - 1)] , 



(92) 



so that 



Z ~ e(f?o) 



N+l 



: exp 



-N 



2A2 



A 



oo 2L 



2(L+fc) 



X [{Qq + 5w) r(L, A;|2L - A;) - OqT{L, k\2L -k - 1)] . 



(93) 



According to classification in seclH ([32|) and ([33]) . the tree and the /-loop approximation 
read 



Ztvec = e ^^°go = exp 



-iV 



(I?o - ^)' 



2A2 



^z-loop = e(^o) 



0^?o)' + A2 



exp 



-N 



{f2o - ujY 
2A2 



A 



l-l ^ 2L 



X [(i7o + Suj) T{L, k\2L - k) - noT{L, k\2L - k - 1)] 



(94) 

2(L+fc) 

(95) 



Up to 3-loop (L < 2), by noting the table[S]in the appendixE]and T{1, 1|0) = -1, T{2, 1|2) 
-3,T(2,2|1) = 35/4,r(2,3|l) = -35/6 for T{L,k\2L -k - 1), we find 



^i-loop = ^(^o) exp 



-N 



(i?o - uY 

2A2 



{^of + A2 



(96) 



^2-loop = ^1-loop + e(^^o) exp 



-N 



(i7o - ujY 
2A2 



N+l 



N 



A^ (I?o - 3(5w) 5A*^ (f2o + Scu) 



(^2o)' + A2 



+ 



(^2o)' + A2 



(97) 



and 



^3-loop = ^2-loop + ^(^o) 



(J?o - g;) 
2A2 



(I?c 



(^?o)' + A2 



iV2 



A^ {f2o - 55w) 



(^?o)' + A2 



7A^ {7f2o + 47(5a;) _ 35A^° {5^0 + DJo;) 385A^2 ^ j^-. 



32 



(^2o)' + A2 



24 



(^^o)' + A2 



72 



(^^o)' + A2 



(98) 
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Since there are two saddle points f?Q ([8^ the total Z is given by 



«p = ^.^op + ^.4<^op' (99) 

where -^^^^QQp has been obtained by putting J7o ^ ■'^o^'* in (195p . 

In the table El we list the result of = 10^(j,6j, IQ-^o; for 10"^ < A < 10 in iV = 2. 
We put (J I— > 1 and write the data of the ratio of Ztroc and -^/.loop^^ ~ ^' '^^ ^'^^ exact 
value. From this, in the weak coupling region, A < 1, the 1-loop approximation almost 
yields the exact value; even in the worst case, = 10"^, only 0.3% error crops up. For a 
whole coupling region including A > 1, the error is within 1.1% under 2-loop and becomes 
venial, < 0.3%, under the 3-loop approximation. 





= 10^ 


A 


Exact 


tree 
(tree) /Exact 


1-loop 
(l-loop) /Exact 


2-loop 
(2-loop)/Exaet 


3-loop 
(3-loop)/Exact 






10~^ 


1.0000 X 10^ 


1.0000 X 10-= 
1.0000 


1.0000 X 10^ 
1.0000 


1.0000 X 10^ 
1.0000 


1.0000 X 10^ 
1.0000 






10~^ 


1.0001 X 10^ 


1.0001 X 10^ 
1.0000 


1.0001 X lO'' 
1.0000 


1.0001 X 10^ 
1.0000 


1.0001 X 10^ 
1.0000 






10-1 


1.0051 X 10^ 


1.0101 X 10^ 


1.0051 X 10^ 


1.0051 X 10^ 


1.0051 X 10^ 






1.0049 


1.0000 


1.0000 


1.0000 






1 


1.5100 X 10^ 


1.8253 X 10^ 


1.5438 X 10^ 


1.5138 X 10^ 


1.5096 X 10^ 








1.2088 


1.0224 


1.0025 


0.9998 






10 


5.2000 X 10^ 


7.5782 X 10^ 


5.3961 X 10^ 


5.1861 X 10^ 


5.1906 X 10^ 






1.4574 


1.0377 


0.9973 


0.9982 










A 


Exact 


tree 
(tree) /Exact 


1-loop 
(l-loop) /Exact 


2-loop 
(2-loop)/Exaet 


3-loop 
(3-loop) /Exact 






10^^ 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 






1.0000 


1.0000 


1.0000 


1.0000 






10-^ 


1.0002 


1.0002 


1.0001 


1.0001 


1.0002 






1.0000 


1.0000 


1.0000 


1.0000 






10-1 


1.0150 


1.0200 


1.0150 


1.0150 


1.0150 






1.0049 


1.0000 


1.0000 


1.0000 








2.5000 


2.8740 


2.4503 


2.4934 


2.4998 






1 


1.1496 


0.9801 


0.9974 


0.9999 






10 


1.5100 X 10^ 


1.8498 X 10^ 
1.2250 


1.4410 X 10'' 
0.9543 


1.5030 X 10^ 
0.9954 


1.5124 X 10^ 
1.0016 










A 


Exact 


tree 
(tree) /Exact 


1-loop 
(1-loop) /Exact 


2-loop 
(2-loop)/Exaet 


3-loop 
(3-loop) /Exact 






10"^ 


1.0001 X 10"^ 


1.0001 X 10^^ 
1.0000 


1.0001 X 10-^ 
1.0000 


1.0001 X 10-^ 
1.0000 


1.0001 X 10-^ 
1.0000 






10"^ 


1.0101 X 10^^ 


1.0101 X 10^^ 
1.0000 


1.0100 X 10-^ 
1.0000 


1.0100 X 10-^ 
1.0000 


1.0100 X 10-^ 
1.0000 






10-1 


2.0050 X 10"^ 


2.0100 X 10^^ 
1.0025 


2.0002 X 10-^ 
0.9976 


2.0049 X 10-^ 
1.0000 


2.0050 X 10-^ 
1.0000 








1.0150 


1.0774 


0.9311 


1.0046 


1.01S2 






1 


1.0615 


0.9173 


0.9897 


1.0002 






10 


1.0051 X 10^ 


1.1105 x" 10^ 


9.1577 X IQI 


9.9945 X lOl 


1.0084 X 10^ 






1.1049 


0.9111 


0.9944 


1.0033 



Table 2: Result of Method(I): = 2, = 1 and cuq = 10^, 1, IQ-^ for 10^^ < A < 10. The 
error is within 1.1(0.3)% under the 2(3)-Ioop approximation for a whole coupling region. 



•Method(II): in ([80]), the saddle points are given by the gap equation, 

A 1 A 



yielding to 



1 



uj + Xtc N ujQ + Xtc 







A2 







ilc — + ^tc 



with 5uj being given by 



This is a cubic equation of f2c, 



(100) 



(101) 



(102) 
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contrary to Method(I), where it was quadratic, (j83p . We write 
The stabihty condition is fulfilled, 

for any (three) saddle points Qc- The values of f{t) and derivatives at Qc are given 



/c = ^ (^^c - t^) - In i7c - + 5^) 



/W = (-r(n-l)! 



AV 1 / A 



n > 3 



(103) 
(104) 

(105) 
(106) 



Then from (|35 



(nj + 3)! 



-Ae(r2c)e(-^c + (^w)]"^ 



n,-|-3 



(f2e + <5c^)"^+' + ^ (r2c)"^+' 



rij + 3 



{Q,f + A2) (I?, + 5uof + f (17,)' 



(n3+3)/2 



(107) 



to give, with using the condition "^1^1=1 ""^j = 2L — 



In exp 



N{Q^-ujy 



2A2 



L=0 fc=0 



(I?,)' + A2) (r?, + 6u;f + § {f2,f 



-iL+fc 



all possible {nj}i-i 



1 



nj+3 



rij + 3 

Z, ([78]) . is expressed, therefore, by a product of a prefactor "P and a loop factor £, 



with 



V = exp 



iv(a-^)' 

2A2 



Z = Vx C 



°° 1 2^ / Nfc 



;^2{L+fc) 



(108) 



(109) 



(110) 



L=0 k=0 

X E n n ■ -I 



(f?,)2 + A2) {n, + <5c^)2 + f (r2,)2 



L+fe 



Ti/ + 3 

all possible {nj}J=i 



1 



(111) 
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Now solve the gap equation to find that there are three kinds of 1/iV series, 



^; = 1,2,3) 



(112) 



where 



a 



; i = l,2 
—5u) ; i = 3 



(113) 



with i?Q^^ being given by (HH), then up to 0{l/N^), 



a 



>(i) 



5,; 



(114) 



Ai = Sn^'' - cu + 6uj ; Bi = 3[^QPj - 2iuj - 6Lo)n^'' - uj6lo - . (115) 

These are sufficient under the 3-loop approximation (0(1/A^^)). Accordingly write Z as 
in (fT09|) such that 



with the prefactor. 



Zii) =p{i) X £W ; {i = 1,2,3) 



where 



fW(i) 



exp 



-iV 



2A2~ 



F»(2) = 



1 + 



+ 



-1/2 



(116) 
(117) 

(118) 
(119) 



and the loop factor, up to 0(1/A^^), 



oo , 2L 



L+k Ej=i"i=2i-A: k 



L=0 fc=0 

1+ 



3^2^(4;i)(5) 5/a2f(3;0(5) 



all possible {nj}^— ^ 

\ 2" 



+- 



AfFW(4) 
(iVFW(4))^ 



+ 



4 FW(4) 6 I F»(4) 



5A2f(6;0(5) A4F(3;i)(5)ir{5;i)(5) 105 /^A2F(4;i)(5)^ 



2 F»(4) 



(FW(4))' 



32 I F»(4) 



105 A6 (F(3;^) (5)) ^ (5) 385 / A^F^S;^ (5) 



(FW(4))^ 



+ 



72 I F{^)(4) 



(120) 
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where 



Note that 



2 / X 2 



M 



; M = 3,4,...,6 . 



1 f2\'^ 
1+ ^ 



+ 



1 



12 



(3) 



1 + 



1 a 



(3) 



; ^ = 1,2 



i = 3 



(121) 
(122) 



(123) 



so that the contribution of the third saddle, i = 3, starts from 0{1/N), that is, 2- loop. 
Therefore we write 



-p(«) 



V, 



1 H 1 — 



(3) 



N 



1 + 



(3) 



N 



i = 1,2 



i = 3 



(124) 



and 



(3) , -^1 



(3) 



(125) 



N 



; i = 3 



ExpHcit forms of those functions, Vq^ ~ 'P^'' and £g ^ ~ d"^^ , are (after lengthy calcula- 
tion) given in the appendix O 

The tree part is, from (|169p and (jl72p . 



7& 



exp 



N 
'2A2 



a 



N 



(126) 



Since ^t^e = 0' ^^i^ equivalent to the one, in Method(I). The 1-loop part reads 



(i) 

1-loop 



' -^0 



e ( I2q*^ ) exp 



N ( I2f ^ - w 



N+l 



,(127) 



and Z| Jq^p = with the aid of (jlSSh and (jl98p . which again matches with (j96p . Next 



(0 

2-loop 



7^, 



1 + 



iV 



1-loop 



1 + 



(128) 
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where 



ii) 



2A2 



1 ( 

2 



(i) 



a 



A2 



f^f) 1 



(0 



+- 



,(0 





SA^ 




SA'^ 








2 + 




3 


4 




6 







from (fTTnll . (fT72]) . (fTTHIl and (fTMl) . For i = 3 we obtain 



(129) 



^(3) 

^2-loop 



(3)^(3) 







e(i7f^) exp 



2-loOp 

TV ( n^^'^ - uj 



1 13 



A^A N12 



(130) 



in view of (fTMl) and (|20T]) . 

Finally Z under 3- loop for i = 1,2 is 



7(i) 

^3-loop 

i-loop 



7^, 



Ar2 



^(i) ^(i) 
_|_ 1 1 j 2 2 11 



N 



(131) 



whose functions 7^i*2'^i,2 are given in the appendix O (fT9T]) . (fTMl) . (fT99]) . and (1200]) . 
Meanwhile z = 3 is 



^(3) ^ ^ 1 
^3-loop iV 



1 + 



(3)- 



(3)- 



^0 +^ 



(3) 



3-loop 



£(3) , ^(3)^(3) 
(3) , -^1 -^0 ^ 2 



iV 



3-loop 



e[nf^ ) exp 



1 313 1 

+ 



iV288 



TV ( Qf^ - UJ 



2A2 

13 
'12 



\ 



eVA^A 



(3) 



+ 



(3)\2' 



1 ; 



A2 



288 2 



(132) 



where 7^2'^^ is given in ()192p and use has been made of ()203p and ()204l) . For Z^^ee ^^^id 
■^i-loop' there is no difference from Method (I). However, in 2- and 3-loop, the third 
saddle starts contributing to give 



^Total ^^(1) ^(2) ^(3) 

/-loop /-loop /-loop /-loop 



(^ = 2,3) , 



(133) 
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whose numerical results, when N = 2 and uj = 1 with ujq = 10^ and 1, are equivalent to 
those of Method(I) for 10"^ < A < 10, which is in the table [2l In the table El we list 
the result of ujq = 10~^, in which we see disparities at 2- and 3- loop in the weak coupling 
region 10~^ < A < 1. Discrepancies are notable, reaching to ~ 600 times to the exact 
value at A = 10^^ in 3-loop. 







A 


Exact 


tree 
(tree) /Exact 


1-loop 
(1-loop) /Exact 


2-loop 
(2-loop)/Exact 


3-loop 
(3-loop) /Exact 






1.0001 X 10"^ 


1.0001 X 10"^ 
1.0000 


1.0001 X 10"^ 
1.0000 


1.0001 X 10-^ 
1.0000 


1.0001 X 10"^ 
1.0000 




10~^ 


1.0101 X 10-^ 


1.0101 X IQ-^ 


1.0100 X 10-^ 


9.6026 X 10-^ 


1.0335 X 10"^ 




1.0000 


1.0000 


0.9507 


1.0232 




10-1 


2.0050 X 10-^ 


2.0100 X 10-^ 
1.0025 


2.0002 X 10-'= 
0.9976 


-3.1834 X lO-'^ 
-15.877 


1.2398 X 10^ 
618.36 




1 


1.0150 


1.0774 


0.9311 


0.8719 


1.2693 






1.0615 


0.9173 


0.8590 


1.2506 




10 


1.0051 X 10^ 


1.1105 X 10^ 
1.1049 


9.1577 X 10^ 
0.9111 


9.9931 X 10^ 
0.9942 


1.0085 X 10^ 
1.0034 





Table 3: Result of Method(II): N = 2, u = 1 and uq = 10'^ for lO'^ < A < 10. Discrepancies 
to Method(I), in the table [21 are acknowledged in 2, 3-loop for 10~^ < A < 10. 



The reason can be seen from ^^^^ (I130p by putting <C 1; 

.(3) 



z 



3-loop 



exp 



2A2 



1 
A ' 



(134) 



which has some peak around A ~ tJo- The graph is shown in the figure [21 implying a large 
deviation at the value u>q = 0.01. Although the case would not be included in the real 



N=2 , (0=1 




Figure 2: The graph of Z 



the value of Z 



(3) 



(3) 

3-loop 



is shown with N = 2,uj = 1 for < A < 1. The vertical line is 



itself. The solid, dotted, and dashed line designate uo = 0.01, 0.05, 0.1. The 



3-loop 

deviation from the exact value (which is almost zero) is seen around A ~ 1.5 when c^o = 0.01. 

physical situation, u- and d-quarks are lighter than s-quark, Method(II) is worse than (I) 
in the above situation. Therefore, a recipe for an approximation Method(I) is better than 
(II) and moreover simpler for actual calculations. 
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5 Discussion 

Under Method(II), we encounter the same situation in the bosonic four-body model under 
the weak couphng region when the "mass", ujq, is tiny: consider cr ([8|) in the introduction 
as bosonic variables, to have 



where in view of ([!]) and 



dy 



ujQ + iXy 



exp [-Nf{y)] 



(135) 



f{y) = ^ + \n{u: + iXy) . 
(We have omitted irrelevant factors, 2tt,N.) Therefore 



dt exp 



-Nf{t) 



fit) = /(t) + - In g{t) ; g{t) = loq + iXt . 



(136) 

(137) 
(138) 



(Again we have switched, y i— > t.) 

Take, for the time being, f{t), fit), g{t) as generic, in other words, start from to 
make a general discussion: there emerge additional saddle points, tQ, in the gap equation 
when N i-^ oo, 



= f'{t) + ^'-^^^git)fit)=0; 
other than /'(to) = such that 

Therefore the additional saddle point (s) is expanded as 

1 



around which 



Now recall that one of the 2-loop terms is given as ( (j35p with ()36| 



2-loop 



(2) 



: exp 



-Nfc 



f< 



?(4) 



(2) 



exp [-Nf, 



By noting 



J(n) ^ An) ^ 



\ 9c 



+ 



{9c 



(139) 



(140) 



(141) 



(142) 



(143) 



(144) 
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we find 



2-loop 



(4) 



?(2) 



3/2 



exp [-iV/, 



{9c? 



exp [-Af/e 



(145) 



around the additional saddle point, which reads in the bosonic case ()138p 



9. 



(1) 



■exp [-iV/cJ ~ , 



(146) 



,(1) 



since gc ^ gc ' ~ A (|138p . again implying a large deviation. 

In summary, Method(II) was superficially simpler than (I) but needs a rather cum- 
bersome procedure in the actual calculation and moreover always seems to suffer from a 
large deviation when in a weak coupling region when the one "mass" ojo is tiny. 

As the final comment, we check the validity of Method (I) in an alternative way: the 
case oi lvq = Lo corresponds to the A'^ = 3 version of the model[8]. 



Zks = / d^'id^'e exp 



(A 



KSj 



2N 



So if we put 



(A 



KSj 



1^' 



(147) 



(148) 



this should agree with our model(A^ = 2). Applying a usual AFM (that is, from the 
expression ([T]) to ([7])), we have results up to 3-loop, which is listed in the tabled! By 
comparing this with the table [2] (of loq = uj), there is no big difference: almost all data 
show that the standard treatment of the table [4] yields a slightly better value except the 
tree approximation in 10"^ < A < 10 where our model under Method(I) results better. 



A 


Exact 


tree 
(tree) /Exaet 


1-loop 
(1-loop) /Exaet 


2-loop 
(2-loop)/Exaet 


3-loop 
(3-loop) /Exact 




1.0000 


1.0000 
1.0000 


1.0000 
1.0000 


1.0000 
1.0000 


1.0000 
1.0000 




1.0002 


1.0002 
1.0001 


1.0002 
1.0000 


1.0002 
1.0000 


1.0002 
1.0000 


X 10-1 


1.0150 


1.0224 
1.0073 


1.0151 
1.0001 


1.0150 
1.0000 


1.0150 
1.0000 


V 2 


2.5000 


3.0574 
1.2230 


2.5433 
1.0173 


2.5032 
1.0013 


2.4996 
0.9998 




1.5100 X 10^ 


2.0188 X 10^ 
1.3369 


1.5484 X 10^ 
1.0254 


1.5106 X 10^ 
1.0004 


1.5089 X 10^ 
0.9993 


Table 4: 


Results of the model (]147[) at = 3, 


uj = 1, and A 
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A The table of T{L, k\2L - k) and T[L, k\2L - k - 1) 
defined by eq.(l52l) 



In this appendix, we list the table of T(L, k\2L — k) / (2(L + A;) — 1)!! instead of r(L, k\2L — k) 
itself; since otherwise the values become very large. The range up to L = 14 is needed for 
the calculation of the Gamma function in SecjS] and to L = 4 in SecHl 



k\L 





1 


2 


3 


4 


5 


6 


7 





1 







-1/6 
47/480 

















1 




-1/4 


-1/8 


-1/10 


-1/12 


-1/14 


-1/16 


2 




1/18 


153/1400 
-493/17280 
77/25920 


3349/30240 
-4049/100800 
190261/29030400 


42131/388080 
-59197/1209600 
1722811/163296000 


605453/5765760 
-161453/2910600 
406957909/27941760000 


655217/6486480 
-81158813/1345344000 
668430857/36324288000 
-39726066467/12070840320000 
4489161401/12070840320000 
-13916027/501645312000 
142181689/101583175680000 


3 






-1/72 


4 






1/1944 


6 








-1/7776 


-503/933120 


-143921/116121600 
44021/522547200 
-1699/503884800 
137/1763596800 


-5689699/2612736000 
74863031/376233984000 
-216793/18811699200 
282133/658409472000 


6 








1/524880 


107/4665600 


7 










-1/2099520 


8 










1/264539520 


9 












-1/1058158080 


-643/63489484800 


-1143257/23702740992000 


10 












1/214277011200 


167/1142810726400 


361657/319987003392000 


11 














-1/857108044800 


-29/1645647446016 


12 














1/254561089305600 


197/1131382619136000 


13 
















- 1/1018244357222400 


14 
















1 /416971064282572800 



k\L 


8 


9 


10 














1 


-1/18 


-1/20 


-1/22 


2 


23763863/245044800 


158899519/1707145440 


11098301/124156032 


3 


-232229821 /3632428800 


-1974182737/29640619008 


-176989210093/2581203905280 


4 


3218616617/146459529216 


131020195003/5187108326400 


1429205465892419/50591596543488000 


6 


-59153229587/13076743680000 


-237186996829/40683202560000 
2187240114496949/2471504555520000000 
-11643095464009/125536739328000000 
23485635969143/3389491961856000000 
-754670748947/2002407559004160000 
4768426956641/315379190543155200000 
-414828299/921562569768960000 
64781047/6516763886223360000 
-7370813/45617347203563520000 


-4136077769339/576345369600000 
27225846275010179/22408307970048000000 
-711993190614367/4943009111040000000 
4800125589203449/388250897448960000000 
-28426615868773/36154580926464000000 
61720076581079809/1639971790824407040000000 
-46647858451427/34060952578660761600000 
14267906296747/374670478365268377600000 


6 


449840404627/747242496000000 


7 


-15727453241 /289700167680000 


8 
9 
10 


570158851513/166867296583680000 
-434386633/2844328919040000 
1125369389/230390642442240000 


11 


-142837/1279948013568000 


12 


2278697/1267148533432320000 


13 


-347/17455617552384000 


-147673007/182469388814254080000 


14 


227/1588461197266944000 


563401/296512756823162880000 


8377793029/640467554738031820800000 


15 


-1/1667884257130291200 


-6271 /400292221711269888000 


-1697149/10674459245633863680000 


16 


1/900657498850357248000 


257/3002191662834524160000 


957911/672490932474933411840000 


17 




-1/3602629995401428992000 


-983 / 108078899862042869760000 


18 




1 /2480410751833883860992000 


41/1049909312945559306240000 


19 






- 1 /9921643007335535443968000 


20 






1 /8483004771271882804592640000 



Table 5: The values of T(L, k\2L - k)/ {2{L + k) - 1)!! for L = ~ 10 in secElas well as secJl 



B The proof of the relation (1491) and of (1541) 

In this appendix, we first prove (j49]) : 



2L 2L 

^T{L, k \ <2L - k) =Y^ T{L, k \2L - k) , (149) 
fe=0 fc=0 

where 

T{L,k\<2L-k) = {-)'{2{L + k)-l)\\ Yl Q,\---Qp[ 

all possible {a^} 
ELi Qc=fc ; (P<fc) 

^ (150) 



(yli + S)"^^ {A2 + 3)"^' ■ ■ ■ {Ap + 3) 



Qi 



23 



k\L 


11 


12 











1 


-1/24 


-1/26 


2 


265842403 /3093554464 


20666950267/249864014400 


3 


-4160074481 /59435616240 


-205489485020713/2894174878795200 


4 


2665668285176389/86090742017280000 


15599534488926830887/466772524452089856000 


5 


-8637934672612781/1011831930869760000 


-40541391139263449/4097919320022528000 


6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 


203058721852862831/128059978750704000000 
-841208011558860343/4033495434608640000000 
405347221415512403/20167477173043200000000 
-198721338986302927/136664315902033920000000 
691811729911733023/8609851901828136960000000 
-202677271543968731/59038984469678653440000000 
3374227932094529/29519492234839326720000000 
-40177324205483/13488137221149661593600000 
79944551761007/1315093379062092005376000000 
-2481097459/2561870218952127283200000 
2893597613/242096735690976028262400000 
-101081/896654576633244549120000 
74363/93537375153331647283200000 
-9637/2381194321760528506552320000 
317/22621346056725020812247040000 
-1/33932019085087531218370560000 


40047282902251643797/20117058480110592000000 
-93845936786285052059/327833545601802240000000 
3993321076456790080543/131133418240720896000000000 
-3058025156141685767/1244621448393523200000000 
21345802689252486829/139397602220074598400000000 
-1543928961441550217/206636445643875287040000000 
19789232975365781149/68666880398580095385600000000 
-103383349078912121/11689718924996373381120000000 
296393081576585261/1367697114224575685591040000000 
-22322956680727/5260373516248368021504000000 
14641217114597/220935687682431456903168000000 
-794890043/968386942763904113049600000 
5917683217/740816011214386646482944000000 
-6685391/111122401682157996972441600000 
56741/165594167212627603331481600000 
-61/42861497791689513117941760000 


22 


1/35272333838948488701496197120000 


347/85508688094420578670293811200000 


23 




-1/141089335355793954805984788480000 


24 




1/175232954511896091869033107292160000 



k\L 


13 


14 



1 








-1/28 


-1/30 


2 


192066102203/2409402996000 


5733412167187/74530866009600 


3 


-316144708749043/4410171243878400 


-5230982263631/72535711248000 


4 


1039787017585599191/29173282778255616000 


86993393978615452781/2310523996037844787200 


5 


-1321164558963178417729/117626676161926643712000 


-409853074642336437191/32674076711646289920000 


6 


124759585247032919932373/51461670820842906624000000 


515724686858368491454901/178997115898584023040000000 
-26661050735382530348701583/55578604486510339153920000000 
15680328374701897384850153/261546374054166301900800000000 
-91936525080313650485417/15932710316247588864000000000 
12636723159765370307618303/28980761501553509007360000000000 
-40157529586233540547571/1529540190359768530944000000000 
23571578581358615769453103/18501318142591760150298624000000000 
-83901448111763933230733/1669291862489482118823936000000000 
884333215764851072621/545074893874116610228224000000000 
-1186820857702705430201/27572773822767445821515366400000000 
1167709777320002702963/1240774822024535061968191488000000000 
-526519789986178363/31019370550613376549204787200000000 
266111014317840767/1054658598720854802672962764800000000 


7 


-3511001126573317919533/9317374453945958400000000 


8 
9 
10 
11 
12 
13 
14 
~l5 
16 


9940709833698231428297/227610147374965555200000000 
-41120299919631976751017/10621806877498392576000000000 
136358632264218467927423/509846730119922843648000000000 
-1026261079863684476899/70256391518917597593600000000 
2015526363139860379999/3161537618351291891712000000000 
-880105003537523330687/39277455587987814560563200000000 
338122116167654792537/530245650437835496567603200000000 
-26876946960034477/1823596152299434247454720000000 
1032930932083247/3730083038794297324339200000000 


17 


-302359614816397/71583162809107792036626432000000 


18 


31701146999281/608456883877416232311324672000000 


19 


-4588191221/8889792134572639757795328000000 


-40513735354423/13142668691752190617924612915200000 


20 


9262374577/2280231682517882097874501632000000 


576706185137543/18728302885746871630542573404160000000 


21 
22 
23 


-7598699/304030891002384279716600217600000 


-2275833187/9120926730071528391498006528000000 


5619821/47884865332875524055364534272000000 


1596747533/985060086847725066281784705024000000 


-13723/33861440485390549153436349235200000 


-17489411/2106934074646523058436039507968000000 


24 


29/29954351198614716558809078169600000 


14256551/436135353451830273096260178149376000000 


25 


-1/700931818047584367476132429168640000 


-211/2213468899097634844661470828953600000 


26 


1/1025112783894592137433843677659136000000 


37/191163223103886645675308844318720000000 


27 




-1/4100451135578368549735374710636544000000 


28 




1/6974867381618804903099872382792761344000000 



Table 6: The values of T{L, k\2L - k) / {2{L + k) - 1)!! for L = 11 ~ 14 in secElas well as secH 

and r(L, k \2L — k^ in RHS is given in terms of the conditional sum '^^=1 Qa^a = 2L — k 
instead of ^^=i Qa^a < 2L — k. When L = (which also means k = 0) the relation 
(fTl9]) holds 

r(o,o I < o) = r(o,o |o) , (151) 

since '^^=1 QaAa < is nothing but X^^^^ QaAa = for any positive Qa,Aa. Therefore 
we can set L 7^ in the following. Likewise when k = 2L 

T{L,2L \<0) =T{L,2L \0) . (152) 

Therefore the target relation ()149p turns out to be 

2L-1 

^T{L,k\<2L-k) = , (153) 

k=0 



24 



with r(L, k\<2L-k) being given by the conditional sum ^ -j Q a.'^a ^ 'ZL — k in 

pSOp . which is further rewritten as 

LHS of (HMD = r(L, I < 2L) + r(L, 1 |2L - 2) 

2L-2 

+ ^ |r(L, A; I < 2L - - 2) + r(L, k + l |2L - /c - 2) | , (154) 
fc=i 

since with the aid of (an obvious relation) 

T(L, k\<2L-k) = T{L, k \2L - k - l) + T{L, k \ < 2L - k - 2) , (155) 
LHS reads 

2L-1 

LHS of (fT53|) = r(L, I < 2L) + ^ {T{L, k \2L - k - l) + T{L, k\<2L-k-2)] 

k=l 

2L-1 

= T{L, I < 2L) + T{L, 1 |2L - 2) + ^ T{L, k \2L - k - l) 

k=2 
2L-1 

+ ^ r(L, A: I < 2L - - 2) , (156) 
fc=i 

whose last term becomes ^^^=2^ ^(-^' ^ I < 2L — /c — 2) , because T(L, 2L — 1 | < — l) = 
0, giving the third term of (jl54p . Meanwhile the third term in (|156p gives the last term 
of (fT54]) by shifting k^k + l. 

The first term cancels the second in RHS of (jl54p : since from ()150p 

T(L,0 I < 2L) = (2L- 1)!! , (157) 

and 

r(L, 1 |2L - 2) = -(2L + 1)!!^^^-!^ = -{2L - 1)!! , (158) 

from pSUp (with changing the conditional sum to QiAi = 2L — 2) and Qi = \ { obtaining 
from 

Ea=i Qa = l; P<1). Therefore if we show 

T{L,K + l\2L- K -2) =-T{L,K \<2L- K -2) ; I < K < 2L - 2 . (159) 

we find that RHS of p54p vanishes, accomplishing the proof of (|149p . 
Now prove (jl59p : LHS reads 

T{L, K + l\2L-K -2) = (-)^+^ (2(L + K) + 1)!! 

Qi!---<3p! {Ai + 3)^5^ ■■■{Ap + 



all possible {Aq} 
Y:1=iQc=k+i ; (P<i^+l) 

In view of (|150p . the conditional sum of RHS, Y.^=i Qa^a < 2L- K -2, is fulfilled 
by putting some Qf^ ; {P £ a) to Qp — 1, which brings the sum Yla=iQa = K to 



25 



Yla=i Qa-'^ = K, while keeping J2a=i QaAa = 2L - K - 2, to give 

T,^=iQcyA^=2L-K-2 

T{L, K \ <2L-K-2) = (-)^ (2(L + K) - 1)11 ^ 

all possible {Aq} 
ELi Qc.=i^+i ; {P<K+i) 
1 1 
^ Qi! • • • {Q(3 -!)!••• Qp! (Ai + 3f' ■■■{Ap + Sf^-^ ■■■{Ap + 3)^^ 

RHS of ()16ip is further rewritten as 



(161) 



RHSof ([l6l]) = (-)-^(2(L + i^)-l)!! ^ 

all possible {Aq} 



13=1 



1 



E^=iQc=A'+l ; {P<K+1) 
1 



(_)^<(2(L + K) + 1)!! 

Er=iQ"^c,=2L-X-2 

E 



1 



1 



Qil---Qpl{Ai + 3f'---{Ap + 3f 



(162) 



all possible {a^} 

where use has been made of 

p 

Y,Qp {Ap + 3) = 2{L + K) + 1 , 



(163) 



/3=1 



obtained from X/^=i ~ — iiT — 2 and X]a=i — 

K + 1. The relation has 
been proved. □ 

Next we show the relation ()54p directly: to this end, prove a general formula 



^p 



dt5(t)e"^^W '^"^ / die-^/W . <7(t) = -fit) + 1 , (164) 
(In our case, /(t) = t - lnt,g{t) = 1/t : B = 0, A = oo.) LHS of (fTM]) becomes 



LHS of (HMD 



p 



whose second term is integrable to yield 



dtf'{t)e 



B 



dfe 



-Nf 



f{B) 



B 



N 



^-Nf{A) _ g-7V/(P) 



(165) 



(166) 



implying no 1/A^ terms for any value of Nf(A), Nf(B). We have proven (jl64p . □ 
In our case, from A = oo, B = 0, that is, from /(oo) = /(O) = oo, 

dt ^ e-^^*^'-^*) , (167) 

which designates r{N + 1) = Nr{N) by multiplying both sides by A^^. 
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>(0 



>(3) 



,(z) 



.(3) 



C Calculation of V'q ' - V\ ' and C'^' - in ( fT24l ) 
and ([1251) 



First expand F(*)(r); (r = 1,2,3), (fTTHI) (fTTOD . such that 

2" 



exp 



iV ( f?^*^ - io 



with 



F(5'^(l)^exp 



2A2 



Ar2 



(168) 



F»(l) 



f2 



2A2 



2 ( - iv) Q'i^ 



(169) 



,F«(1): 

■ (2) Fj'^(2) ' 

^ iV ^ Af2 



F«(l) 



with 



F«(2)^ ff2( 



Fo«(2) 



F«(2)^ 



2A2 

TV 



Ff(2) 



^^Q*^") exp 



(2) 



1 



1 

















1^ . (170) 



(171) 

(172) 
(173) 



(3) 

(Note that there is no need for F2 (r); {r = 1, 2).) 



F«(3) = 
with 



1 + 



ii) 



+ 



-1/2 



F«(3) 



N iV2 



F«(3) 



Ff)(3) 



iVA 

A2 



^2l 



a 



a 



(3) 



(3) 



2A2 
27 



1 
2 

1 / 12 



a 



(3y 



, (174) 



(175) 



(176) 



F«(3)^ 



3(Ff^(3) 



+ 



(i) 



+ 



'0 



(Again no need for fP{3).) Further F(*)(4),([l2l]), for i = l,2 reads 



F»(4) =F(S^)(4) 



with 



F»(4) 



Meanwhile for i = 3 



212 



- + r^^r^ + 



J7, 



+ A2 



7VF(3)(4) = Fo^^^(4) 



i7(3)(4) 



with 



(4)^2^ + 



(3) 



/2 



(3)- 



(3) 



a 



(3) 



+ 



(3) 



i?(A-/;i)(5);(M = 3,...,6), (fT22D . for i = 1,2 reads 



A2 



with 



M 



(5) ^ (^f?^'^ + 5a; j , (5) ^ 

Then for i = 3 

with 



+ - 



Ma 



^(A^;3)(5y 



M = 3,... ,6 , 



i.(M;3)(5)^^^(3) 



Af 



n(3) 



(3) ■ 



28 



Therfore the prefactor reads 



exp 



N ( - w 



N+i e(f?, 



for i = 1,2 and the leading term of i = 3 is 

21 



V{ = exp 



~2>? 



(188) 



(3) 



'1 



where use has been made of the relation 



^ — + ^ 



A2 



eVNX 



i = l,2 
-1 i = 3 



(189) 



(190) 



in view of (I113p and the gap equation (I83p for i = 1, 2 or ()114p with ()115h for i = 3. Then 



r=l 



(3), 



i7, 



(3) 



r=l 



a 



(3) 



(191) 

(192) 



and 



+ E i^i^^^(r-')i^f^(r) + -^^E^i^^W- (193) 



(0 



r=l 



I7g + Jli; r'>r=l 



i7Q + (5w r=l 



Now calculate terms in the loop factor (jl20p up to the 3-loop approximation: 
A2 



= 1 + 



+- 



iVF»(4) 
(iVF«(4))^ 



.^i?(4;.) + ^f^(3;i)^' 

4 6 V 



5A 



lOSA^ 



_R(6;i) + 7A4i?(3;i)ij(5;i) + 
2 32 



72 



where we have introduced 



FW(4) ° 
whose expansion coefficients are 



1 + 



N 



{M;i) _ -t^o 



(M;i) 



(5) 







FoW(4) 



+ 



M-2 



1 

iV2 



(194) 



(195) 



ni;^) +A2 

M-2 



(196) 



a 



(3) 



A2 



29 



R 



{M;i) 



(i) 



(M-2)-i 



(3) 



J7 



(3)- 
1 



(3) 



(3) 









2 






A2 





^2^^+ 5lo) 



I?f] + A2 



With the aid of these, 2 = 1,2 part is given by 



C 



(i) 




1 



3X^ 



/^f) + A2 



+ 



3A_2 
4 

5A6 



{4;i) 



+ 



5A4 



(3;i) 



I^f)+A2 



n 3 



.(197) 



(198) 



(199) 



'-2 



of) + A2 



5A2 



?(6;i) 



i4 Tj(3;i) Tj(5;«) 



105A^ 



(4;i)^ 



105A6 







2 385A8 / 



(3;i)\" „(4;i) 



^0 + 72 



(3;i) 



3A^ E!f^I^^ 5A6(^?fV^ 



+■ 



385 



A6 



f2«) + A2 
Ai2 



n 3 



+ 



329 



A8 



4 F(^)(4) 
105 



+ 



3 

Aio 



32 



A2 



I??)+A2 



•A 2 



H 5 



72 



+ A2 



£W^W(4) _ /i + ^1A_L^ . (200) 



4 Fo»(4) 



Fo«(4) 



For i = 3, in view of (jlSip . a naive 1/A^ expansion is broken down; all factors in 
front of the loop expansion, (jl20p . behave 0{1). Nevertheless we take up to the terms 
of 0{l/ {NF^^\A)]) under the 2-loop approximation, and of o{l/ {NF^'^\A)Y^ under 
3-loop so that 



(3) 



1 + 



2-loop 



(4) 



13 
12 



(201) 



30 



with the use of pl2]) and and for 0(l/iV) terms we obt 

A2 



iam 



,(3) 



-loop ^o'^ (4) 



-— ^0 



,(4;3) ^ |^^{3;3)^ ^ 



6 



Ff(4)) + 



7?(4;3) „{4;3) 



?(3;3) 



(3)' 



+ 



A2 



(202) 



where use has been made of (fT82]) . (fT83]l . ([TM]) . and (fTW]) . 

In a similar manner, we obtain the 0(1) term in 3-loop, such that 



,(3) 



c 



(3) 



3-loop 



+ 



loop(F(^^(4) 



5A2 (6;3) 7x4p(3;3) „(5;3) 105A^ / (4;3) 



105A^ 



, , p(3;3) 



313 

288 ' 



(203) 



and for 0{1/N) terms 



(3) 



A^ 



3-loop 2-loop (Fo^^H4) 

105A*^ 



„(3;3)\2 (4;3) , 385A8 / (3; 



-2Ff^(4) 



+ |_^46;3)^(6;3) ^ 7;,4^(3;3)^(5;3) ^^(3;3) ^ ^(5;3)^J ^ 1^ 



16 



.(4;3) 



R 



.{4;3) 



105A^ 



.(3;3)V p(4;3) 



■0 



2^(3;3)^^(4;3)N 



.(4;3)\ , 385A*^ 



18 



(3;3) 



R 



{3;3) 



. .X 137 (3;3) 117 (4;3) „„(5;3) 5 (6;3) 



13 
12 



a 



(3) 



=0 



A2 



(204) 
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